Classiness and constrained dynamics of a short-range non-disordered spin model 
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We study the low temperature dynamics of a two dimensional short-range spin system with uniform 
ferromagnetic interactions, which displays glassiness at low temperatures despite the absence of 
disorder or frustration. The model has a dual description in terms of free defects subject to dynamical 
constraints, and is an explicit realization of the "hierarchically constrained dynamics" scenario for 
glassy systems. We give a number of exact results for the statics of the model, and study in detail 
the dynamical behaviour of one-time and two-time quantities. We also consider the role played by 
the configurational entropy, which can be computed exactly, in the relation between fluctuations 
and response. 



I. INTRODUCTION 

Understanding the nature of the low-temperature dy- 
namics of glasses and other strongly interacting many 
body systems remains one of the outstanding operiprob- 
lems in condensed matter and statistical physics . 

Consider, for example, the case of glass-forming super- 
cooled liquids. One of the reasons for glassiness in such 
systems is that the rearrangements of atoms necessary for 
their relaxation involves activation over energy barriers. 
It is natural to plot the logarithm of the viscosity or equi- 
libration time against the inverse temperature, giving a 
so-called Arrhenius plot, which should take a straight-line 
form if the barrier heights remain constant with varying 
temperature. For "strong" hquids such as Si02, the 
Arrhenius plot is indeed a straight line, but for "fragile" 
liquids it is not, following instead the empirical Vogel- 
Fulcher law exp[const./(T — Tq)], although other forms 
not displaying a finite T singularity can be fit too j^]. If 
the dynamics of fragile liquids is due to activation (which 
is not the only possible explanation) , this implies that the 
energy barriers grow with decreasing temperature, pre- 
sumably because they increase with the increasing size of 
correlated regions in the system. 

Another set of interesting questions, most relevant for 
fragile liquids, revolves around the possible existence of 
an ideal (continuous) phase transition to a true ther- 
modynamic glass state at some temperature Tk (the 
Kauzmann temperature) lying below the glass temper- 
ature UJ^. Although the existence of such a transition 
would resolve Kauzmann's paradox, in which the extrap- 
olation of the configurational entropy of the supercooled 
liquid appears to pass below that of the crystal j|] , and is 
also supported by the analogy between fragile glasses and 
discontinuous mean field spin-glasses 0, there is some 
evidence against it. For example, it has been shown nu- 
merically that a thermodynamic phase transition is 
absent for a system of polydisperse hard-disks, a typical 
model glass former. It has also been argued that the ex- 
trapolations which yield a positive Tk in fragile liquids 
composed of molecules of finite size are fiawed |9| . 

Given the many interesting open questions in this field. 



it is important to find simple and if possible solvable 
models for the various features displayed by glasses. The 
discontinuous mean-field spin-glasses, such as the spher- 
ical p-spin model pO[ | , provide a good example of such a 
model system. However, while some dynamical features 
of these models [| 11 are also observed in simulations of 
Lennard- Jones glasses [|2|,^ , the presence of disorder in 
the models and, crucially, of long-range interactions lim- 
its their usefulness as models of systems which are intrin- 
sically short-ranged and disorder-free, like hard spheres 
in two and three dimensions. 

An alternative approach to modeling these latter sys- 
tems is the "hierarchically constrained dynamics" of 
Palmer et al. |lj], in which it is hypothesized that, for a 
strongly interacting system displaying glassy behaviour, 
it should be possible to describe the dynamics in terms 
of hierarchies of degrees of freedom, from fast to slow, in- 
dependent of the presence of disorder or even frustration. 
These hierarchies would be weakly interacting in the en- 
ergetic sense, but their dynamics would be constrained, 
the faster modes constraining the slower ones. It is known 
that the presence of kinetic constraints in the dynamics 
can directly induce glassiness. A good example is the 
kinetically constrained Ising chain p5|-p7[|, for which the 
Hamiltonian is trivial, but the transition rates for the 
flipping of individual Ising spins depend on the states 
of neighbouring spins. Kinetic constraints can also arise 
as a result of entropic barriers, as in the Backgammon 
model, for instance [Q. Fragile glass behaviour has also 
been observed in kinetic lattice-gas models [^,^. Al- 
though interesting in their own right, these models have 
somewhat ad hoc dynamics, and represent only one side 
of the scenario discussed in Ref . | . A recent compara- 
tive study of one-dimensional constrained kinetic models 
has been given by Crisanti et al. pl| ]. 

In this paper, we study in detail a model introduced 
by Newman and Moore ||2^, which is a form of two- 
dimensional Ising model with uniform short-range fer- 
romagnetic interactions. Despite the absence of either 
disorder or frustration, this model displays glassy be- 
haviour at low temperatures. The cause of this behaviour 
is the presence of energy barriers which grow logarithmi- 
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cally with the size of correlated regions. The model has 
a dual description in terms either of strongly interact- 
ing spins subject to simple single-spin-flip dynamics, or 
of free "defects" subject to a constrained dynamics. It 
is thus an explicit realization of the constrained dynam- 
ics scenario of Palmer et al. [0. Other non-disordered 
short-range spin systems displaying glassy features are 
three-dimensional Ising models with competing nearest 
and next-nearest neighbour interactions [ p3[ | , or with fer- 
romagnetic four-spin plaquette interactions p^-|2^ . 

The paper is organized as follows. In Section || we 
describe the model and extend the exact solution of the 
statics given in Ref. |Q by calculating all equilibrium 
spin correlation functions, and showing that the phase- 
transition to an ordered state occurs only at zero tem- 
perature. In Section [II, we give our principal results, 
which concern the out-of-equilibrium dynamics of the 
system following a quench from a random configuration 
to low temperatures. We show that the equilibration 
time diverges with an exponential inverse temperature 
squared law, similar to that found for the asymmetri- 
cally constrained Ising chain. We study the behaviour of 
one-time quantities, like energy density and spin corre- 
lations and two-time quantities, like autocorrelation and 
response functions. We also study the out-of-equilibrium 
fluctuation-dissipation relations and relate them to the 
configurational entropy, which can be calculated exactly. 
Our conclusions are given in Section IV. 



II. MODEL AND STATIC SOLUTION 

We consider the model introduced in Ref. |^2|, which 
consists of Ising spins a = ±1 on a triangular lattice 
with uniform short-range three-spin ferromagnetic in- 
teractions: each spin interacts only with its nearest- 
neighbours, and only in groups of three lying at the 
vertices of a downward-pointing triangle on the lattice. 
Note that this is distinct from the model of Baxter and 
Wu p^ ], which has interactions on upward-pointing tri- 
angles also (and which is not glassy). The Hamiltonian 
for the model is 



(1) 



where the indices m and n run along the unit vectors of 
the lattice ai = x and 0,2 = \{x + VSy)- The constant 
in the Hamiltonian is added for convenience to make the 
minimum possible energy equal to zero. The model can 
also be formulated using the defect variables 



in terms of which the Hamiltonian is 



H 



(2) 



(3) 



On lattices of size a power of two in at least one direction, 
with periodic boundary conditions, there is a one-to-one 
correspondence between spin and defect configurations, 
and hence the partition function is given by 



Z = {2e-i'^ coshi/3) 



(4) 



where N is the number of spins, and we set J = I from 
here on. The equilibrium energy density is then 



(1 + (t» = i(l -tanhi/3). 



(5) 



We now wish to invert Eq. (H) and write the spins as 
functions of the defects. It is known |^ that if the spins 
are represented by the variables Smn = ^{o'mn + 1) G 
{0, 1} the spins below an isolated defect on the lattice 
form a Pascal's triangle mod 2 (i.e., a triangular array 
with the binomial coefhcients ( " ) mod 2 as entries) . This 
means that the value of the spin is given by the 
superposition mod 2 of the Pascal's triangles of all defects 
Tki — +1 with m ^ I < k < m and / > n. The inverse of 
the transformation (^ then reads 



O'rnn — 



(6) 



-L<k<r, 



With this result we can compute the equilibrium corre- 
lation functions of the spins. First, the magnetization is 
given by averaging (^): 



{cr„in) = -(tanhi/3) 



where 



E (™-lO-od2 



(7) 



(8) 



is the total number of ones in the inverted Pascal's tri- 
angle with its tip at site (m, n). This number diverges 
faster than the linear size of the triangle, so in the ther- 
modynamic limit we obtain 



-1 




for T = 
for T > 0, 



(9) 



for all (m, n), implying that the system has a T = static 
phase transition. 

Arbitrary correlation functions can be computed in a 
similar manner: 



■Grn'ri') = -(tanhi/3) 



(10) 



where now Mmn---m'n' is the total number of ones in 
the superposition mod 2 of the inverted Pascal's trian- 
gles starting in positions (m, n) through (m',n'). Since 

(0) ~ (;) ~ ^'^^^ ^""^"^ right edges of a Pascal's 

triangle contain only ones, so any superposition of two 
triangles always has an infinite number of ones coming 
from the edges. This implies that all 2-spin correlations 
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vanish at T > 0. The first non-zero correlations are those 
for three spins at the vertices of inverted equilateral tri- 
angles of side 2*^. In this case A/" = 3'' , and the correlation 
is 



exp 



1 



2T2 ln2 



(14) 



C^. — (Cmn (T„ „_|_2'! Cr„j_2'= ,ri+2'« ) — ^ (tanhi/3)^ 

Notice that 

C^^i(tanh i/3) = Cf ^(tanh^ i/3). 



(11) 



(12) 



A similar scaling relation is seen in the one-dimensional 
Ising model. 

At low temperatures the system is in a scaling region, 
and using Eq. (11) we find a correlation length of 



e= (lncothi/3)-'"2/in3^ 



(13) 



III. DYNAMICS 

Newman and Moore found that the model stud- 
ied here shows glassy behaviour under a single-spin-flip 
dynamics. They showed that following a quench from 
T ~ oo the system was unable to equilibrate in finite 
time at low enough temperature (T < 0.2 for tmax ~ 10^ 
in their simulations). The system also fell out of equi- 
librium for exponential cooling with a variety of cooling 
rates. The reason for this glassiness is to be found in the 
details of the system's dynamics. Single-spin fiips cor- 
respond to flips of the defect variables on triples of sites 
forming upward-pointing triangles on the lattice, and sets 
of such flips can be combined to flip upward-pointing tri- 
angles of side 2'"' for any integer k. Any isolated such 
triangle with fc > is locally stable; in order to remove it 
we have to cross an energy barrier of height k. Thus as 
the system relaxes it has to cross energy barriers which 
grow logarithmically with the size of the equilibrated re- 
gions. 

Another way of looking at this is to observe that at low 
temperatures the flipping of an isolated defect is heavily 
suppressed, since such a flip requires the creation of two 
other defects and so incurs a net energy penalty. Thus, 
a defect requires the presence of another neighbouring 
defect to be able to flip, a situation highly reminiscent 
of the facilitated kinetic Ising model of Ref. ||l^ ; the de- 
fects are non-interacting in the Hamiltonian, but their 
low temperature dynamics is effectively constrained. 

At low temperatures excitations of linear size larger 
than 1 can only be annihilated via activation. From the 
observation that barriers grow logarithmically with lin- 
ear size it is straightforward to estimate equilibration 
time. The rate of relaxation of an excitation of lin- 
ear size d is given by the Arrhenius formula T{d) ~ 
exp(— lnd/rin2). Thus after time t the average linear 
distance between defects is d ~ i^'"^. The equilibrium 
value of this distance at low T is deq ~ exp(/3/2) [see 
Eq. (^] , and hence the equilibration time is 



This exponential inverse temperature squared (EITS) 
form is similar to that obtained in Ref. p^ ] for the asym- 
metrically constrained Ising chain (ACIC) except 
for the 2 in the denominator which is due to the two- 
dimensional nature of our model. The EITS form has no 
finite-temperature singularity, which is consistent with 
the fact that our model has no finite-temperature phase 
transition. The behaviour of the relaxation, however, is 
"fragile". In fact, the EITS form can be approximated 
by a Vogel-Fulcher form with Tq ~ ^Tg. A system which 
has an analogous but "strong" behaviour is the Frenkel- 
Kontorova model at low temperatures above its 

T — phase transition, the relaxation of defects is hin- 
dered by the presence of constant energy barriers, and 
the equilibration time diverges in an Arrhenius way. 

Equation (|lj) clarifies the observations made in 
Ref. [Q. In the case of relaxation following a quench 
from Tq = cx) to T, the highest temperature at which the 
system was found to equilibrate within time <max = 10^ 
was T = 0.2, while for T — 0.18 or smaller the system 
was unable to reach equilibrium. This is explained by the 
fact that ioq(T = 0.2) - 10*, while toq(T = 0.18) - 10^''. 
In the case of exponential cooling, our expression for icq 
implies a maximum allowable cooling rate of 7max(2^) ~ 
l/tcq(r) if we wish to remain in equilibrium at tempera- 
tures T and greater. For the different annealing simula- 
tions of ||2^ , the temperatures at which the energy ceases 
to follow the equilibrium curve are given by the solution 
of 7max(?^) — 7 for the cooling rate 7 used. 



A. One time quantities 

We now consider the behaviour of one-time quantities 
for the model, such as energy density and spin correla- 
tions, following a quench from T = cx) to a low tempera- 
ture. 

The T — 00 configuration is random, and therefore has 
zero correlation length. As the system relaxes following 
the quench it tries to increase its correlation length to the 
equilibrium value, Eq. (p^, which is finite but large at 
low temperature. In terms of the defects, the system is 
trying to decrease the defect density or equivalently the 
internal energy per spin e from its initial value of Eq = i . 
In Fig. we show numerical results for the evolution in 
time of the energy density after a quench from T = 00 
to a variety of final temperatures. The simulations were 
performed on a 256 x 256 rhombic lattice using a Bortz- 
Kalos-Lebowitz continuous time algorithm ||29| . As the 
figure shows, after an initial temperature-independent 
exponential relaxation corresponding to the (barrierless) 
removal of pairs of neighbouring defects, the energy dis- 
plays "plateaus" , which are more pronounced the lower 
the temperature. These plateaus are the result of the 
system becoming trapped in locally stable configurations. 
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FIG. 1. Energy density e = {H)/N as a function of time 
following a quench from T = oo to T = 0.10, 0.12, 0.15 and 
0.20, from Monte Carlo simulations for a system of 256 x 256 
spins. Time is in Monte Carlo steps per spin, and all magni- 
tudes are in dimensionless units. Inset: the same data as a 
function of rescaled time i/ = Tint. 

As we can see, the time taken by the system to escape 
these plateaus becomes larger with decreasing tempera- 
ture. 

As shown in the previous section, the typical length- 
scale in the system (i.e., the typical distance separat- 
ing defects) grows with time as t'^ , and in the absence 
of other length-scales we might thus expect the data in 
Fig. |l| approximately to collapse when plotted against 

, or equivalently against v — Tint on the logarithmic 
scales used in the figure. (A collapse of this kind was 
found for the ACIC in Ref. 0.) In the inset of Fig. 
we show that this is indeed the case for our model. Save 
for the initial exponential transient, the collapse of the 
curves for different temperatures is very good. Moreover, 
the first plateau seems to extend from = to 1, the 
second from = 1 to 2, and so on (a behaviour also seen 
in the ACIC). 

At low temperatures, the relaxation of the system can 
also be regarded as an anomalous coarsening process. For 
the case considered here of periodic boundary conditions, 
there exists a unique minimum energy configuration of 
the model in which all the spins point down and there 
are no defects. For free boundary conditions, however, 
there are 2^^^^ degenerate minima. (The spins along 
the bottom and right hand side of our rhombic system, 
for instance, may be chosen arbitrarily, with the rest be- 
ing uniquely fixed by the requirement that there be no 
defects.) Following a quench, spins start to rearrange 
locally to eliminate defects, regardless of the boundary 
conditions, and thus form domains of the various free- 
boundary minima. The imperfect matching of these do- 
mains will be marked by the presence of defects. 

The coarsening dynamics for T close to zero can 
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FIG. 2. Average distance between defects d as a function 
of rescaled time v — T\n{t). The dot-dashed line corresponds 
to the T = coarsening approximation. 

be studied approximately using the generating function 
method of SoUich and Evans [|7| . Consider a description 
of a configuration of the system in terms of domains con- 
taining no defects (r = —1) bounded by sites containing 
defects {t — 1). In contrast with the one-dimensional 
model studied in , such a description cannot be made 
exact for our model. Nevertheless, as far as average prop- 
erties go, such as distance between defects, we can think 
of low energy configurations as an approximate tiling of 
defect-free parallelograms each delimited by a defect in, 
say, its top right-hand corner. This effectively maps the 
problem into one dimension. As we now show, this crude 
approximation, which allows us to apply the method of 
Ref. [0, gives reasonable results. 

When T — > 0, we have Scq and the plateaus of 
Fig. |l| become distinct stages in the dynamics. Dur- 
ing stage n, all domains of linear size d in the interval 
2''~^ < d < 2^ are annealed away. The time-scale for 
this process is 0(e'^'^), and differs by a factor of from 
the time-scale associated with the following stage. In 
this limit, the master equation for the coarsening pro- 
cess can be solved using generating functions. If Pk{d) 
is the probability distribution of lengths d at the be- 
ginning of fc-th stage, Gfc(z) = Y.2''-^<dPk{d)z'^ is the 
probability generating function for this distribution, and 
Hk{z) = X]2'=-i<(i<2'= Pk{d)z'^ is the equivalent function 
for the active domains — those which will be annealed 
away during the current stage — the following recursive 
relations are obtained [|l^ 

Gk+i(z) = 1 + [Gk{z) - l]exp(i/fe(z)). (15) 

Given an initial distribution Pa{d), these equations can 
be solved numerically. 

In Fig. ^ we show the average linear length of the do- 
mains d{t) — 1/ \/s(t) as a function of i/ = T In t from 
our simulations. As T ^ 0, we expect d{t) to become 
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FIG. 3. Two-spin correlations for various times as a func- 
tion of distance. 



a "staircase" function, with the steps extending between 
integer values of ly. The approximate heights of these 
steps are given by the derivative of the generating func- 
tion at z = 1: {d)k = J2d'^^k{d) = G'j,(l). These steps 
are also plotted in the figure (dot-dashed line) and, as we 
can see, the agreement is quite good. 

We now turn to the behaviour of the spin correlation 
functions as we approach equilibrium. In equilibrium, as 
shown above, the magnetization and all two-spin correla- 
tions vanish at any finite temperature. Moreover, at all 
times following a quench, the average magnetization and 
all equal-time 2-spin correlations vanish as well, as shown 
in Fig. ^. This is a consequence of the three-spin interac- 
tions: regardless of the value of a spin i, the Hamiltonian 
favours a neighbouring spin j equally to be up or down, 
since the only interaction between the two spins also in- 
cludes a third spin of unknown orientation. 

In equilibrium, the first non-zero correlations are those 
for triplets of spins at the vertices of downward pointing 
triangles of linear size 2*^. This is also the case in the out- 
of-equilibrium regime following a quench. In Fig. W we 
show the absolute value of the 3-spin correlations | | 
as a function of time after a quench, for various values 
of k. The behaviour of these curves illustrates the nature 
of the stages by which relaxation takes place. Initially, 
correlations at all length-scales are zero. The system has 
no barriers to relaxation of excitations of length-scale 1 , 
so jCp I starts to grow immediately after the quench. At 
J/ = Tlni = 0, \C'[^^\ grows exponentially fast to a first 

(3) 

quasi-stationary value, while all C^>]^ remain zero; this 
first plateau corresponds to partial equilibration up to 
length-scales of 2. At = 1, the system starts to relax 
up to length-scales of 4, with C2 becoming nonzero. 



while C, 



(3) 
k>2 



0. And so forth. 



Within the plateaus all one-time quantities have 
roughly stationary values, and so it is natural to ask to 



(3) 

FIG. 4. Three-spin correlation functions as a function 
of time, for linear sizes 2'"^'''i'^'^'*. We plot the absolute value 
of the correlations. 



what extent the system is in a quasi-equilibrium state. 
One answer to this question comes from examining the 
agreement between our non-equilibrium correlation func- 
tions and the exact relations (^2|) for the equilibrium 
correlation function. In Fig. ^ we compare |Cq'^^| with 

I^F^I'"^^'^ and find that the out-of-equilibrium corre- 
lation functions still collapse approximately when scaled 
appropriately, the collapse getting better as time pro- 
gresses. This may indicate that the later stages can be 
described by an approximate equilibrium at an appropri- 
ate effecti ve tem perature. We discuss this point further 
in Section III C. 



B. Two-time quantities 

We now turn to the behaviour of two-time quantities 
in the out-of-equilibrium regime following a quench from 
high temperature. We first consider the local two-time 
spin autocorrelation function 



C(i, t„ 



N 



{t)(ymn{tw))- 



(16) 



In Fig. 1^ we show C(i, iu,) for quenches from T — 00 
to T = 0.12 for three different values of the weighting 
time tw, as a function of the rescaled time difference i> — 
TlnT, with T = t — t.^. The three values of used 
correspond to = 0, 1, 2, i.e., to the starting poi nts o f 
the first three of the plateaus discussed in Section III A . 
In Fig. 1^ we present the same correlations as a function 
of t/t^. 

Since the system does not reach equilibrium at T = 
0.12 on the timescales simulated, we do not in general 
expect the correlation functions to be functions of the 
time difference t = t — only. In fact, as Fig. shows. 
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FIG. 5. Three-spin correlations scaled according to 
Eq. 

the behaviour of the correlation functions has a clear de- 
pendence on the waiting time. On the other hand, as we 
can see from Fig. neither does it obey a simple aging 
form, scaling with t/t^, as observed in other cases 
This is because the length-scale associated with the re- 
laxation of the spins does not grow as a simple power 
of time (see Figs. || and ^). In general, if equilibra- 
tion is associated with the growth of a length-scale l{t), 
then the two-time correlation functions should scale with 
l{t)/l{tw) Behaviour of this kind was found in the 
autocorrelation functions for the ACIC j2^, where the 
appropriate length-scale is the average distance between 
upward-pointing spins. In our model, however, rescaling 
the spin correlation functions by the distance between 
defects fails, and we have not been able to find a suitable 
length-scale such that the scaling with l{t)/l{tw) holds. 

Another observable of interest is the response func- 
tion tu})^ which measures the response of the spins at 
time t to a small field applied at time . Comparison of 
such a response function with the two-time correlations 
studied above can reveal violations of the fluctuation- 
dissipation theorem (FDT) which are expected in sys- 
tems displaying aging. In general, in order to obtain 
the response function corresponding to the local correla- 
tion function calculated above, we would need to apply a 
field to only a single spin on the lattice, or equivalently we 
could apply a random field and measure the staggered re- 
sponse [|o| . Both of these approaches present significant 
numerical challenges, and require a substantial invest- 
ment of CPU time in simulation to extract clean results. 
For our model, however, this turns out t o be unnecessary 
because, as mentioned in Section [II A, all off-diagonal 
spin-spin correlations vanish, implying that two-time au- 
tocorrelations of the magnetization m = N^^ l^mn^'mn 
are proportional to the local spin correlations thus: 



FIG. 6. Local spin-spin correlation functions as a func- 
tion of the scaled time difference Tlnr, following a quench 
to T = 0.12. Inset: the same for T = 0.20. 




FIG. 7. Local spin-spin correlation functions as a function 
of t/t^, for T = 0.12. Inset: the same for T = 0.20. 

This means that we need only to measure the response 
to a uniform field. 

We have performed simulations in which the Hamilto- 
nian was perturbed with a small uniform magnetic field 
applied at a time following a quench to finite temper- 
ature: AH{t) — hod{t — tiu) J2mn '^mn- Wc measured the 
integrated response (or relaxation function) 



x(i,i« 



1 



y^(0'nm(^))fe, 



(18) 



(17) 



which is equivalent to measuring the change in magneti- 
zation per spin of the system, divided by the strength of 
the applied field. We have experimented with a variety of 
different values for the field strength h^; the results pre- 
sented here are for ho = 0.02, which we find to be well 
within the linear regime. In Fig. ||we show the measured 
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FIG. 8. Integrated response function vs. rescaled time dif- 
ference, for T = 0.12 and TXnt^ = and 1. Inset: the same 
for T = 0.20. 

response for T = 0.12 and Vyj = Tlnt^, =0 and 1. No- 
tice that the response is not in general a monotonically 
increasing function of time, but instead shows cusps at 
times corresponding to the jumps between plateaus. For 
higher temperatures the cusps are less pronounced and 
eventually disappear (see inset). 

A similar non-monotone behaviour is seen in the re- 
sponse function for the defects in our model, and can be 
understood physically by the following argument. The 
energy density e plays a role for the defects equivalent to 
that played for the spins by the magnetization to, and the 
equivalent of a small applied magnetic field 5h is a small 
change in temperature 5T — 6h. Consider then the 
behaviour of the energy density shown in Fig. |l]. At 
low temperatures we expect e to follow a staircase func- 
tion, with steps extending between integer values of v. 
When the field is applied (i.e., when the temperature is 
changed), the perturbed energy density is another stair- 
case function, with each step having a slightly higher 
(lower) value for negative (positive) field, and lasting a 
longer (shorter) time, the steps now extending between 
integer values of i^(l — TSh). The corresponding inte- 
grated response is given by the difference between the 
perturbed and unperturbed energy densities, divided by 
the field. Within each plateau this difference is mostly 
small and roughly constant, but a sharp change is seen at 
the boundaries between steps. The unperturbed {T — * 0) 
density jumps to the next plateau when v is an integer, 
but the perturbed one only does so when ^{l — TSh) is 
an integer, and thus the absolute value of the response is 
large during the short interval between these two events. 
Once the perturbed density jumps to the next plateau, 
the response becomes small again. This increase and de- 
crease at integer v is responsible for the humps seen in 
the integrated response. The argument generalizes to 
perturbations with a staggered random field. 



The response function for the ACIC was studied in 
Ref. ||2^, where it was found to be a monotonically in- 
creasing function of time. However, as mentioned above, 
the defect representation of our model behaves very simi- 
larly to the ACIC, and we therefore conjecture that non- 
monotone behaviour will be seen in the ACIC also at suf- 
ficient low temperatures, for precisely the reasons given 
above. 



C. Conflgurational entropy and 
fluctuation— dissipation relations 

A configuration of our model is a local energy minimum 
if and only if no two defects occupy adjacent sites on the 
lattice. In the context of glasses, such states are known 
as inherent structures |Q . A nice feature of the model 
is that the set of inherent structures is isomorphic to the 
set of allowed states of Baxter's hard-hexagon model [ p2| , 
allowing the distribution of inherent structures to be cal- 
culated exactly. 

The grand-partition function of the hard-hexagon 
model on a lattice of N sites is 

N 

Zn{z)^Y. 9{M,N)7J', (19) 

where z is the fugacity and g{M, N) is the number of ways 
of placing M non-adjacent particles on the triangular lat- 
tice. In our case, M = Ne corresponds to the total en- 
ergy of our model, so 5jv(e) = g{Ne, N) gives the density 
of states for inherent structures with energy e. For large 
N, g{n, N) is exponential in A^, so that (?Ar(e) = e-'^^''^^\ 
where S'c(e) is the configurational entropy of the model, 
i.e., the entropy density of metastable states. In this case, 
Eq. ( p^ ) becomes 

Zn{z)^ J de exp[N{e In z + Sc{e))]. (20) 

For large N the integral is dominated by the saddle point 
in the exponent, and we obtain 



Sc[e{z)] = In K{z)~e{z) In z, (22) 

where k(z) = limN^cx3[ZN{z)Y^'^ is the partition func- 
tion per site of the hard-hexagon model, which is known 
exactly in the thermodynamic limit [ p2| . Between them, 
Eqs. ( pT| ) and ( |2^ ) determine Sc{£) parametrically. 
At low defect densities, 5'c(e) reduces to 

Sc{e)^ -elne + 0{e), (23) 

which is the general form for a low concentration of non- 
interacting point defects in an ordered structure ||^ . The 
configurational entropies of the disordered Ising chain, 
constrained Ising chains, and the Backgammon model 
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all have this asymptotic form Note that Sc has 

infinite slope with respect to e at £ = where it van- 
ishes, which is consistent with the fact that none of these 
models has a finite temperature phase transition. How- 
ever, since S'^ increases logarithmically as e is decreased, 
an extrapolation of S'^(0) from the behaviour at finite e 
would wrongly suggest a finite slope and therefore a fi- 
nite Kauzmann temperature. It has been argued that 
this mechanism would rule out an ideal-glass phase tran- 
sition in materials composed of limited size molecules and 
conventional molecular interactions 

Having calculated correlation and response functions 
for the spins in our model (Section IIIB), we can use our 
results to study the relation between fluctuations and re- 
sponses, which in some other systems is related to the 
configurational entropy. The natural way to do this is 
by means of a parametric plot of x{tw + T,tw) versus 
C{tw+T, tw) for fixed tw pl . For a system in equilibrium 
at temperature T, such a plot would be linear with slope 
— 1/T, in accordance with the fluctuation-dissipation 
theorem. In a glassy system on the other hand, the 
FDT is normally violated in the out-of-equilibrium low- 
temperature regime as t^, — s- oo, and this violation en- 
codes important information about the system's dynam- 
ical behaviour (see Ref. ^ for a review). The classic ex- 
ample is the mean- field p-spin spin-glass [p^ , for which 
the FDT plot is piecewise linear: for large values of C, 
corresponding to relaxation of fast degrees of freedom, 
the plot has a slope of —1/T as it would in equilibrium, 
but for smaller C, corresponding to the relaxation of slow 
modes, it has slope — 1/Teff, where Toff > T is interpreted 
as a effective temperature for these modes . Further- 
more, in the p-spin model T^g is numerically equal to the 
inverse of the slope of the configurational entropy at the 
asymptotic energy density p5[ |. Similar behaviour has 
also been observed in more realistic model glass form- 
ers ll|Jl|. 

At sufficiently low temperatures our model is clearly 
out of equilibrium, and although it does not reach the 
long-time asymptotic regime corresponding to approach 
to equilibrium in our simulations, it is possible that 
the FDT plots can still provide information on the re- 
laxation process, and maybe that there exist effective 
temperatures associated with the different time-scales in 
the problem |34|. In Fig. ^ we show the FDT plot for 
the spin response and correlation functions. The un- 
usual non-monotonic shape is a consequence of the non- 
monotonicity of the response function. 

The curves should be "read" from right to left in the 
plot, and are composed of a sequence of segments, each 
associated with one of the plateaus. The starting point 
of each segment and the part of the curve in which the 
response is increasing correspond to the time the sys- 
tem spends within the relevant plateau. The maximum 
and the downward portion of the segment correspond 
to the transition to the next plateau. The first part 
of each segment has a shape similar to that found in 
Refs. pGl and p^ for the one-dimensional Ising model 




FIG. 9. Parametric plot of integrated response vs. the 
two-time autocorrelation function of the spins for T — 0.12 
and Tlntw — and 1. The dotted line corresponds to slope 
— l/T. In the inset we show the same for T = 0.20 and 
TlnU = 0. 



following a quench into the temperature scaling region: 
the FDT curve there has slope —1/T at C = 1, and 
bends continuously with decreasing C to reach a slope 
of —1/2T at C = 0. (Note that this behaviour is very 
different from that seen in domain growth in higher di- 
mensions [pO|.) In our model, each of the segments of 
the FDT plot starts with slope approximately —1/T (in- 
dicated by the dotted line in the figure), but the sub- 
sequent shape of the curve varies from one segment to 
another. The similarity within the plateaus to the be- 
haviour of the one-dimensional Ising model is consistent 
with what was found for the statics of the model in Sec- 
tion |l|. 

An important open question is whether the configura- 
tional entropy of Eqs. (21) and ( [2^ ) plays any role in the 
out-of-equilibrium dynamics. A plausible explanation for 
the shape of the FDT curves of Fig. IqI is that the ini- 



tial part of each segment corresponds to quasi-stationary 
thermal excitations of fast modes, while the latter part 
corresponds to slower large-scale rearrangements arising 
from occasional jumps between local minima before the 
transition to the next plateau takes place. If the slope 
of the latter part of each segment corresponds to an ef- 
fective temperature for that segment, we might expect 
these temperatures to be related to the rate of change 
of Sc{s) at the energy density of the associated plateau. 
It is noteworthy that the final slopes of each segment in 
Fig. H are roughly equal to 1/S'^{ek), with values approx- 
imately i, ^, and i for the first three plateaus. This 
observation is however speculative and furthermore con- 
trasts with the conclusions of Ref. for the ACIC, so 
it deserves further investigation. 
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IV. CONCLUSIONS 



activated processes. 



In this paper we have studied the low temperature be- 
haviour of the glassy two-dimensional spin model with 
uniform ferromagnetic short-range three-spin interac- 
tions introduced by Newman and Moore. The model has 
a dual description in terms either of interacting spins or 
of free defects, the mapping between the spin and defect 
representations being one-to-one. This allows us to com- 
pute exactly all equilibrium correlation functions, both 
for spins and defects. We have also shown that the model 
displays no static phase transition at finite temperature. 

Despite the simplicity of its statics, the model's low- 
temperature single-spin-flip dynamics is highly non- 
trivial. A spin flip in the model corresponds to the flip- 
ping of three neighbouring defects, which implies that 
at low temperatures the dynamics of the defects is con- 
strained: a defect can only be flipped if another neigh- 
bouring defect is present. This in turn implies that the 
relaxation of isolated defects is an activated process, and 
the size of the corresponding energy barriers are found 
to grow logarithmically with the distance between de- 
fects. The dual representation of the model in terms of 
strongly interacting spins with a simple dynamics, and of 
free defects subject to kinetic constraints, is an explicit 
realization of the hierarchical constrained dynamics sce- 
nario of Palmer et al. . 

We have studied in detail the dynamics of the 
model after a quench to low temperatures. The pres- 
ence of logarithmically growing barriers leads to non- 
Arrhenius relaxation, the equilibration time being of 
the exponential inverse temperature squared form tcq ^ 
exp(l/[2r^ln2]). One-time quantities such as internal 
energy display "plateaus" in their equilibration profiles, 
which correspond to the trapping of the system in lo- 
cal energy minima. Each plateau is associated with a 
specific stage in the dynamics, the fc"^ plateau corre- 
sponding to the partial equilibration of length-scales up 
to 2*^. The behaviour of observables related to the de- 
fects is strikingly similar to that seen in the asymmetri- 
cally constrained Ising chain. For example, in the T — > 
limit the average distance between defects can be well 
approximated using the analytic methods of SoUich and 
Evans which yield exact results for the asymmetri- 
cally constrained model. 

We have also studied two-time quantities for the 
model, like spin autocorrelations and response functions. 
At low temperatures, the response functions have the 
unusual property of being non-monotonic: they display 
humps at exactly those times at which the system jumps 
between plateaus. This behaviour has also been ob- 
served in other models at times well within the acti- 
vated regime, such as the constrained Ising chains and 
the Backgammon model in one dimension pl| , models 
for two-dimensional froths js^, and vibrated granular 
media ]39| ]. An important open question is whether this 
is a generic feature of the out-of-equilibrium dynamics of 
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